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Abstract 

The time-dependent transmission coefficient for the generalized Kramers problem with ex- 
ponential memory friction has recently been calculated by Kohen and Tannor [D. Kohen and 
D. J. Tannor, J. Chem. Phys. 103, 6013 (1995)] using a procedure based on the method of 
reactive flux and the phase space distribution function. Their analysis is restricted to the high 
friction regime or diffusion-limited regime. We recently developed a complementary theory for 
the low- friction energy-diffusion-limited regime in the Markovian limit [Sancho et al. , to appear 



in J. Chem. Phys.; cond-mat /9806001 1 . Here we generalize our method to the case of an expo- 
nential dissipative memory kernel. We test our results, as well as those of Kohen and Tannor, 
against numerical simulations. 
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1 Introduction 

The classic work of Kramers [jlj] on reaction rates in which the effect of the solvent was taken 
into account in the form of a Markovian dissipation and Gaussian delta-correlated fluctuations has 
spawned an enormous and important literature that continues to flourish ^ ||, [|, |5| ) . The enormous 
literature on the theoretical front since the appearance of Kramers' seminal paper has evolved in 
many directions that include more formal derivations of Kramers' own results, extensions to larger 
parameter regimes and to non-Markovian dissipation models, generalizations to more complex 
potentials and to many degrees of freedom, analysis of quantum effects, and application to specific 
experimental systems. 

One recent direction, developed by Kohen and Tannor (KT) ||, deals with the derivation of 
the rate coefficient in the Kramers problem and in the more general Grote-Hynes problem M (that 
is, the Kramers problem extended to a non-Markovian dissipative memory kernel) so as to obtain 
not only the asymptotic rate constant but the behavior of the rate coefficient at all times. This 
derivation is based on the reactive flux method ]2|, ||, 0], thus paralleling closely and usefully the 
methods used in numerical simulations of the problem. KT analyze in detail the time it takes 
the rate coefficient to reach its stationary (equilibrium) value and the way in which this value is 
approached. Their extensive analysis, however, does not cover a number of parameter regimes, 
nor do they check their time-dependent results against numerical simulations. In this paper we 
complement their work by extending the parameter regime of analysis and checking their results 
as well as our new ones against numerical simulations. 

The non-Markovian generalization of the Kramers problem is based on the dynamical equations 
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where q{t) is the time-dependent reaction coordinate, a dot denotes a time derivative, T(t — t') is the 
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dissipative memory kernel, V(q) is the potential energy, and f(t) represents Gaussian fluctuations 
that satisfy the fluctuation-dissipation relation 



(f(t)f(t')) = k B TT(t-t'). (2) 

k B is Boltzmann's constant and T is the temperature. The potential V{q) is a double-well potential 
that is often (and here as well) taken to be of the form 



vw-*,[*i-\ + i)-?V-V- (3) 

The parameter Vo can be used as the unit of energy, and henceforth we set it equal to unity. The 
barrier height is assumed to be large compared to the temperature (i.e., k B T -C 1/4). 

The problem of interest is the time-dependent rate coefficient k(t) for an ensemble of particles 
whose positions evolve as realizations of q(t). The coefficient k{t) measures the mean rate of passage 
of the ensemble across the barrier at q = from one well to the other. The asymptotic value of 
this crossing rate, k = k(oo), is associated with the rate constant of the process represented by the 
reaction coordinate. One usually focuses on the corrections to the rate obtained from transition 
state theory (TST) and therefore writes 

k(t) = n(t)k TST (4) 

where k TST is the rate obtained from transition state theory for activated crossing, which for our 
problem and in our units is [Q, ||, [U| 

k TST = V2 e -l/4fc fl T_ (5) 

The deviations from this rate constant are then contained in the transmission coefficient «(t). The 
construction of n(t) is discussed in subsequent sections. 

In this paper, as do KT, we deal with the exponential memory kernel 

r(t) = ^ e -i*i/ T (6) 

r 

where 7 is the dissipation parameter and r is the correlation time. The Kramers problem cor- 
responds to the limit r — > 0. In their work, KT derived an expression for K,(t) in terms of the 
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Figure 1: Transmission coefficient K s t vs dissipation parameter 7 for ksT = 0.025 and various 
values of the memory correlation time r. These results are obtained from direct simulation of 
Eq. d]) as described in Section 2.2 , As is well known, the turnover moves to higher values of the 



friction parameter with increasing r because the effective dissipation at given 7 is weaker as r 
increases. Also, the maximum transmission coefficient decreases with increasing r. 

time-dependent phase space density, assumed a particular form for this density, and studied how 
n(t) goes to its steady state value as the phase space density relaxes to equilibrium. In carrying out 
this program, however, they relied on approximations that are appropriate only for the diffusion 
limited regime (large 7) and therefore were able to calculate n{t) only in this regime. 

In a recent paper [10] we developed a theory for the time-dependent transmission coefficient in 



the energy-diffusion-limited regime (small 7) for the original Kramers model. Thus, to supplement 
the work of KT and our own work in the Markovian regime, here we formulate the complemen- 
tary theory for the time-dependent transmission coefficient in the non-Markovian energy-diffusion- 
limited regime and compare our results with numerical simulations in this regime. As in the 
Markovian case [pX| ], our theory captures the complex oscillatory and stationary state behavior of 
n(t) in this regime very well. We also carry out simulations of the time-dependent transmission 
coefficient in the diffusion-limited regime to assess and confirm the validity of the KT formula- 
tion. We will show examples of the good agreement of their results with simulations, both in the 
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"non-adiabatic" regime (large 7 but 7 < t) and in the "caging" regime (large 7 with 7 > r). 

In Section ^ we briefly review the reactive flux formalism and describe the numerical methods 
used in our simulations to implement the formalism. Typical simulation results in various regimes 
are shown in this section. In Section ^ we recall the theoretical methods of analytic calculation 
of the transmission coefficient in the diffusion-limited regime (that is, KT's method) and energy- 
diffusion-limited regime (our method). The generalization of our previous Markovian results to the 
non-Markovian regime is detailed in this section, as are comparisons of these various results with 
numerical simulations. We conclude in Section ||. 

2 Reactive Flux Formalism and Simulations 

2.1 Reactive Flux Formalism 

Two decades ago saw the development of the reactive flux formalism for the rate constant over a 
barrier [^, |8|, H. This formalism was important because it made possible the efficient numerical 
simulation of rate constants without having to wait the inordinately long time that it takes a 
particle at the bottom of one well to climb up to the top of the barrier. In the reactive flux 
method the problem is formulated in terms of the particles of the distribution that are sufficiently 
energetic to be above the barrier at the outset. In this way, it is not necessary to wait for low-energy 
particles (the vast majority of all the particles) to first acquire sufficiently high energies via thermal 
fluctuations. 

We follow the notation of KT. The rate constant k(t) in the reactive flux formalism is 

_ (ep(q )e P [q(q ,v ,t]) _ (v Q 5(q Q )e P [q(q Q ,v Q ,t)]) 
[) (0 R (qo)) (0 R (9i)) ' [) 

where the top of the barrier is at position qo = 0, 9r(x) = 1 if x < and otherwise, and 
9p = 1 — Or. The brackets ( ) represent an average over initial equilibrium conditions. With 
the particular choices made in Eq. (|7]) one is calculating the transition rate from the left well to 
the right. KT proceed through a series of steps that finally yield for the transmission coefficient 
introduced in Eq. (Q) the relation 

K(t) = I ^ f J dv v e- m ^ 2k * T X (v , t) (8) 
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where 



oc 



x (v ,t) = dq dvW(q,v,t;q Q = 0,v ) (9) 

JO J-oo 

and W(q,v,t; qo = 0,vq) is the conditional phase space distribution function that corresponds to 
an ensemble of particles starting at (qo = 0, vq) at t = 0. 

The (nonequilibrium) conditional probability distribution W clearly lies at the crux of the 
calculation: it is the distribution associated with the generalized Langevin equation ([!]). In Section |3| 
we discuss the approximate solution of this problem (which can not be solved exactly) to obtain 
analytic results. The above definitions are sufficient for the numerical simulations. 

2.2 Numerical Simulation Method 

The numerical solution of Eqs. (|l]) is performed according to the following main steps. First, we 



write the problem in the entirely equivalent form[ll, 12] 



P ■ dV(q) . p z 

z, z = - 7 +rj(t) (10) 



m ' dq t t 

where 

(v(tHt')) = 2 -^5(t-t>). (11) 

The integration of this set of equations is carried out using the second order Heun's algorithm [13|, 
which has been tested in different stochastic problems with very reliable results |l4j| . A very small 
time step is used, ranging from 0.001 to 0.0001, as in Ref. ||. The numerical evaluation of the 
transmission coefficient K(t) follows the description of Refs. ||, |9|. We start the simulation with N 
particles (between 1000 and 4000 depending on the circumstances), all of them above the barrier at 
q = 0, half with a positive velocity distributed according to the Boltzmann distribution in energy, 
which translates to the velocity distribution v exp(— v 2 /2/cbT), and the other half with the same 
distribution but with negative velocities. The transmission coefficient is extracted from these sets 
of simulated data by calculating ^ 

, . N+(t) N-(t) . s 

" (() = NW) - ACM' <12> 
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where iV+(i) and iV_(i) are the particles that started with positive velocities and negative velocities 
respectively and at time t are in or over the right hand well (i.e. the particles for which q(t) > 0). 

As discussed by Straub et al. |J, whereas the exact transmission coefficient reaches a constant 
non-zero value that corresponds to the equilibrium transition rate for the problem, the transmission 
coefficient calculated using the reactive flux method flattens out but continues to decrease with time. 
If the temperature is low, this decrease is slow and the results appear flat for rather long simulation 
times. In this case the value of n(t) in this flat region is identified with K s t- If the temperature 
is not so low, then the decaying tail is extrapolated back to its intersection with the vertical axis 
according to the relation n(t) ~ n s te~ Kt where K is the decay constant of the tail. The value of 
the intersection is then identified as K s t- 

2.3 Simulation Results 
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Figure 2: Numerical simulation results for the transmission coefficient k(V) vs time for memory 
correlation time r = 3.0 and temperature h^T = 0.025, and four values of the dissipation parameter 
7. Diffusion-limited non-adiabatic regime: 7 = 2; diffusion-limited caging regime: 7 = 10; energy- 
diffusion-limited regime: 7 = 0.01; intermediate regime: 7 = 0.3. 

Figure ^ exhibits examples of the temporal behavior of the transmission coefficient, Eq. 
associated with the system (|l|), for four typical different parameter regimes. The input parameters 




7 



of the problem are the dissipation parameter 7, the correlation time r, and the temperature ksT. 
The memory correlation time in the figure is held fixed at the value r = 3.0 and the temperature 
at fcgT = 0.025. The 7 = 2.0 curve is typical of the diffusion-limited "non-adiabatic" regime. 
This is beyond the turnover value for r = 3.0 in Fig. |l| but still small enough to be in the non- 
adiabatic regime. The transmission coefficient decays monotonically to the equilibrium plateau. 
The curve differs from the typical behavior in the Markovian diffusion-limited regime [10, |(| in 
that the initial decay is Gaussian rather than exponential || (this detail is just barely visible in 
the figure). The 7 = 10.0 curve is typical of the regime well beyond the turnover dissipation and 
well into the "caging regime" || . The transmission coefficient decays in an oscillatory fashion with 
an oscillation frequency that can be associated with an effective potential discussed later. The 
7 = 0.01 curve is typical of very low dissipation, well within the energy-diffusion-limited regime. 
The transmission coefficient remains at its initial value up to a time beyond which it decays in an 
oscillatory manner to its equilibrium value. As shown later, this behavior is quantitatively captured 
by our theory. The oscillation frequency in this case is completely different from that of the caging 
regime. Finally, 7 = 0.3. is below the turnover in Fig. |] but not much below, and the associated 
transmission coefficient exhibits a mixture of features characteristic of diffusion-limited behavior 
(the initial decay to an apparent plateau) and energy-diffusion-limited behavior (the subsequent 
oscillatory decay to equilibrium). 

3 Theory 

The conditional probability distribution W(q, v, t; qo = 0, vq) associated with Eq. (|l|) can not be 
calculated exactly for the bistable potential (|3|). Various approximations associated with different 
parameter regimes and the resulting time-dependent transmission coefficients are discussed in this 
section. 

3.1 Diffusion-Limited Regime 

KT focus on the regime of moderate to high dissipation, the regime first discussed in detail by 
Grote and Hynes pi] . This is the regime to the right of the turnover in Fig. |lj; the specific values of 
7 included in this regime clearly depend on the correlation time r. To obtain an explicit expression 



for the probability distribution W, KT adapt the theory of Adelman [15] to the assumption that 



the barrier is parabolic and the wells are infinitely deep. With this approximation, KT derive 
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expressions for K(t) from which the Grote-Hynes equilibrium results for the transmission coefficient 
are recovered in the long-time limit. For arbitrary times they obtain the expression 

K (t) = C ^ (13) 

y/cd(t) + 

where 

A n (t) = -^(C 2 v (t)-C 2 (t) + l), (14) 
m 



C q (t) = l+ I dt' C v (t% (15) 
J o 

and 

C v (t) = C- l [s 2 - 1 + st(s)]- 1 = Y, c ^ lt - ( 16 ) 

i 

Here £ _1 denotes the inverse Laplace transform, T(s) is the Laplace transform of the memory 
kernel T(t), and the /ii are the roots of s 2 — 1 + sT(s) = 0. For the exponential memory kernel (||) 
r(s) = 7/(1 + rs); the equation is then cubic so there are three roots ^j. The largest one of them 
is positive and it determines the equilibrium transmission coefficient. The others may be real or 
complex, depending on the parameter values, and this in turn determines whether the transmission 
coefficient approaches equilibrium in an oscillatory or monotonic fashion. KT show that in the 
Markovian regime (r — > 0) the decay of n(t) is monotonic and exponential. On the other hand, 
for r large compared to the time scale ^1 of the reaction (the precise limits on r for this to be 
the case are discussed in KT), they show that in the so-called "non-adiabatic" regime the decay 
is monotonic, whereas in the "caging" regime the decay is oscillatory. This behavior is easy to 
understand: for a time interval smaller than r one can roughly approximate T(t — t') as constant. 
Performing the integral in Eq. (|l|) over such a time interval, one obtains an effective potential V e ff 
of the form 

Veff(q) = V(q) + l£ (17) 
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If 7 < r (non-adiabatic regime), this is still a bistable potential, albeit with a smaller barrier than 
that of the "bare" potential V(q). However, if 7 > r (caging regime), then the effective potential 
over this time interval is monostable with a frequency oo ca gi n g ~ (7/T- l) 1 / 2 about q = 0. 

In Fig. H we show again two of the simulation results shown in Fig. [2], those for which the KT 
theory is valid, as well as the results of Eq. (|b|). The monotonic "non-adiabatic" decay and the 
oscillatory "caging" decay are clearly captured very well by the theory, including the oscillation 
frequency in the latter (which for these parameters is uj cag i ng = 1.53) and the equilibrium values in 
both cases. 




Figure 3: Numerical simulation results and theoretical results for the transmission coefficient n(t) 
vs time for memory correlation time r = 3.0 and temperature k^T = 0.025, and two values of 
the dissipation parameter 7 in the diffusion-limited regime. Solid curves: simulation results for 
7 = 2.0 (diffusion-limited non-adiabatic regime) and 7 = 10.0 (diffusion-limited caging regime). 
The dashed curves in each case are calculated from the theory of Kohen and Tannor as given in 



Eq. ©. 



The KT theory is not appropriate for low dissipation parameters and indeed it fails in that 
regime because it does not take into account the inertial recrossing of the barrier that dominates 
this regime. Thus, a different approach is required in that case 
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3.2 Energy-Diffusion-Limited Regime 

When the dissipation is very slow, that is, well to the left of the turnover in Fig. |], the dynamics 
of the system is dominated by the slow variation of the energy. Thus, a particle that starts above 
the barrier may recross the barrier many times before it loses sufficient energy to be trapped in one 
well or the other. We note that there are interesting questions concerning the role and applicability 
of the reactive flux method and hence of the transmission coefficient in the prediction of time- 
dependent reaction rates in the low dissipation regime Here we do not deal with these issues. 
Rather, we take the reactive flux method as our starting point and carry out our analysis within 
that framework. 

Because of the slow energy variation it is here convenient to rewrite the dynamical equations 
(FJ) in terms of the displacement and the energy (henceforth we set the mass m equal to unity) 

E =Y + v ^ ^ 



instead of the momentum [16, |17|. This simple change of variables leads to 



q=p(t) (19) 



E = -p(t) f dt' T(t - t') p(t') + f(t)p(t) (20) 
J o 

where p(t) is understood as a shorthand notation for {2[E(t) — V(q(t))]} 1 ^ 2 . This set can of course 
not be solved exactly either, but one can take advantage of the fact that for slow dissipation the 
temporal variation of the energy is much slower than that of the displacement. 

When dissipation is very slow, the particles that start above the barrier lose their energy very 
slowly as they orbit around at an almost constant energy during each orbit (recrossing the barrier 
many times if they start with sufficiently high energy). One can calculate the approximate energy 
loss per orbit and thus keep track of how long it takes a particle of a given initial energy above 
the barrier to lose enough energy to be trapped in one well or the other. There is a distribution of 
such particles above the barrier, and the time that it takes the ensemble to lose enough energy to 
be trapped in a well is correspondingly distributed. 

The principal ingredients in the calculation of the transmission coefficient are: 1) The time t e 
that it takes an orbiting particle of energy e above the barrier to complete a half orbit, that is, to 
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return to q = 0; 2) The energy loss in a half orbit, yu(e), of a particle of initial energy e above the 
barrier; and 3) The distribution of times at which particles eventually become trapped in one well 
or the other. 

To calculate the time that it takes a particle of energy e above the barrier {E = 1/4 + e) to 
complete one passage from q = to the edge q+ of the potential well and back to q = 0, we assume 
that the energy remains fixed during this passage. g+ is given by the solution of E = V(q), i.e., by 
q + = (1 + \J4lE) 1 / 2 . Independently of the value of r and 7, the time t e for this half-orbit is then 
simply obtained by integrating Eq. (|i~9|): 



t = 2 r dq (2D 

- Jo {2[E-V(q)]} 1/2 

With the potential (||) this can be expressed in terms of a complete elliptic integral of the first kind. 
For low temperatures the particles that determine the behavior of the transmission coefficient are 
primarily those just above the top of the barrier. For these energies an excellent approximation to 
the complete elliptic integral yields ||, 1C] 



16 , 16 

t £ wln — =ln— . (22) 



The corrections to Eq. (22) are of O(elne). 

The energy of the particle of course does not in fact remain constant as the particle orbits 
around. The calculation of n{e) will be postponed for the moment - we return to it later in this 
section. The result does of course depend on 7, t, and temperature. Suffice it for now to say that 
to an acceptable level of approximation the energy loss fi can be taken to be independent of e. 
Although this assumption is not essential, it does simplify the calculations considerably. 

The steps in the construction of the transmission coefficient from these components have been 
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detailed in our earlier work[|l0|1. Here we simply display and interpret the final result: 



K (t) = 1 + 2 f2(-l) n e~ Mt)/kBT , < t < In H 

n=l 



/i 



l_2e-M/*BT + 2 g ( _ 1) n e -/n(i)/fc fl r ln^<t<21n- ' (> 

n=2 



/i (2!)V2 M > 



= i _ 2e -^ T + 2e~ 2 ^ T + 2 f](-l)« e ~/«(')/^ T , 21n— ^L- < t < 31n^y^-, 

(23) 

and so on. The function £ — fm(t) is the solution of the relation t — t £ -\- t e —^ + • • • + te—(m~i)^,i 
which upon exponentiation turns into the m th order polynomial equation 

(16) m e -* = £ [e - n] ■ ■ ■ [e - (m - (24) 

The solution of Eq. (^) can in general not be found in closed form for m > 3. However, an excellent 
approximation is 

f m (t) » \m - (miyH p. + 16e-'/ m (25) 

This form is exact for m = 1, and it is exactly correct for all m at the upper limit of the time range 
that defines the trapping of the particles whose initial energy is e = mix. In other words, f m (t) as 
given in Eq. (^5|) is exactly correct at the particular time t = tn + t2u + • • • + t mfl . The asymptotic 
limit of Eq. ( |2"3| ) gives for the equilibrium transmission coefficient 

K st = K(t -> OO) = 1 + 2 ^(_l)« e -»MA fl T = tanh / \ (26) 

n=l ^ ' 



It is useful to examine the information contained in the various contributions to Eq. (p3[). The 
first line contains all first recrossings of the barrier by all particles that start with sufficient energy 
above the barrier to recross is at least once [e > fi). The function f\{t) accounts for the fact 
that this first recrossing occurs at different times for particles of different initial energy, the last 
ones (those closest to the barrier that do make it around) recrossing at time t = ln(16//i). The 
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particles of higher energy recross first because they are orbiting at higher velocity, and may be back 
for their second recrossing while those of lower energy are still awaiting their first crossing. The 
time distribution of the second recrossing is contained in f2(t)- Even while this is occurring, some 
particles might already be undergoing their third recrossing, as contained in fs(t). There are fewer 
and fewer of these faster particles because of their initial thermal distribution - this information, 
too, is contained in the exponential factors. Of course while all these recrossings are going on, the 
particles are losing energy and, depending on their initial energy and how often they have gone 
around, they become trapped at sequential times as expressed in the subsequent lines of Eq. (p3|). 
In constructing these results we have averaged over a half orbit in calculating the time for the 
return of a particle to the origin, and have assumed that in the calculation of this period the energy 
loss of the particle over the half orbit is negligible. Thus, we have approximated the evolution of 
the energy of a particle as determined by Eqs. (|i~9| ) and (p0|) with f(t) = by simply assuming an 
energy loss of fj, at the end of each half orbit. 

The only remaining calculation is that of the energy loss jx per half orbit. This quantity must 
be calculated from the generalized Langevin equation Eq. (f20|). 

3.2.1 Markovian limit 



In our earlier work [ 10 1 we detailed and discussed the calculation of // for the Markovian case. We 
showed that when the fluctuations are ignored entirely, the equation for the rate of change of the 
energy, E = — 27 [22 — V(<z)], averaged over a half-orbit leads to the approximate equation 

47 1 

E . (27) 

£.5-1/4 



Integration then directly leads to an implicit algebraic equation for fi(e) |T0|. It turns out that fi 
does depend on the initial energy e of the particle above the barrier - specifically, fi increases with 
increasing e, but that this dependence is mild, as seen in Fig. |] for two values of the dissipation 
parameter. An upper bound is the simple relation between the energy loss per half orbit and the 
dissipation parameter 

V = \l- (28) 

This result agrees with that obtained via weak-collision arguments or small dissipation arguments 
by a variety of essentially equivalent routes. The total variation of fi with e is of order 10% for 
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dissipation parameters in the energy-diffusion limited regime, as can be seen in the figure. 




Figure 4: Energy loss per half orbit as a function of energy in the Markovian limit (r — ► 0) for 
two values of the dissipation parameter: 7 = 0.005 (dashed curve) and 7 = 0.01 (solid curve). 

Thermal effects tend to reduce /x since thermal fluctuations allow the particle to gain some 
energy as it orbits. We refer the reader to Ref. [jl0| for a discussion of thermal effects. 

Figure ^ shows the transmission coefficient n{t) as a function of time in the Markovian limit 
(r — ► 0) for a typical set of parameters in the energy-diffusion-limited regime. The distinctive 
features of the time dependence are 1) the time at which the transmission coefficient drops rather 
abruptly from its initial value of unity, and the slope of this drop; 2) the frequencies and amplitudes 
of the oscillations; and 3) the asymptotic value, identified as the equilibrium value n s t of the 



transmission coefficient. The dashed curve is the result of our theory, Eq. (23), with the bare value 
/i = 47/3 for the energy loss per half orbit. We stress that there are no adjustable parameters in 
these curves. The agreement between the theory and simulations is clearly very good, and typical 
of parameters in this regime. Without adjustable parameters, the theory captures each of the three 
distinctive features listed above. Indeed, we stress that with a single expression we are able to 
reproduce the temporal behavior of the system over essentially all time scales. 

The theoretical equilibrium value of the transmission coefficient achieved by the dashed curve 
in Fig. H is a little higher than the simulation result. We also show the equilibrium transmission 
coefficient obtained when we include thermal effects [|H]] (thick short line that intersects the right 
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Figure 5: Transmission coefficient of an ensemble of particles at temperature ksT = 0.025 and 
dissipation parameter 7 = 0.005 in the Markovian limit. Solid curve: simulation of 4000 particles. 
Dashed curve: our theory, Eq. (p3|), with the bare value = 47/3 = 0.006667. Thick line that 
intersects the right vertical axis: value of the equilibrium transmission coefficient K s t obtained from 



Eq. (26) using the value of renormalized by thermal fluctuations, which for this temperature is 
H(T) = 0.004907. 



vertical axis). This value of the stationary transmission coefficient agrees extremely with the 



simulation result. In [10] we discuss the connection of our equilibrium result to other predictions 
S 



3.2.2 Non- Markovian limit 

To calculate /x in this regime we again begin by ignoring the fluctuations and thus approximate the 
energy evolution by the first portion of Eq. (|20|), 



E = -p{t) I dt' T(t - t') p{t'), (29) 
Jo 

where p(t) = {2[E(t) - V(q(t))}} 1/2 . To integrate this expression over a half-orbit is complicated 
by the fact that the two p factors appear with different time arguments. One requires an analytic 
expression for pit) in order to carry out this integral. It turns out that an excellent approximation 
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to p(t) is obtained with just two Fourier components of the form 

p(e, t) A(e) cos 

J + B(e) cos • ( 30 ) 

We have explicitly indicated the energy dependence in p and in the coefficients. The latter are fixed 
by imposing two conditions. One is that the initial value of p be the correct one associated with 
its definition and the fact that q(0) = 0, 

p(e, 0) = A + B = V2e. (31) 

The other is the equal-time average over a half orbit, already calculated for the Markovian case [cf. 
Eq. (f27f: 

1 f ts 9,^ A 2 + B 2 4 

d„ ^ = — = 5C- < 32 > 

One obtains 




Figure ^ is a typical example of the agreement between the exact deterministic pit) obtained via 
numerical integration and the approximation just presented. In this figure we have chosen an energy 
e associated with our subsequent discussion of the transmission coefficient. 

With this analytic expression for p(t) one can now explicitly calculate the energy loss [i per half 
orbit. Substitution of Eq. ([}(]) into Eq. fl29|) , and integration over t between and t e with fixed e 
immediately leads to the expression 




with A and B given earlier. This result is shown in Fig. |j] for r = 3.0 and 7 = 0.01. 

The differences in Figs. || and [7| warrant discussion because of the "opposite" behaviors of 
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Figure 6: Momentum as a function of time for a particle moving in the potential (|3|) with energy 
e = 0.003175 above the barrier. Solid curve: exact result. Dashed curve: approximation (|30|) with 



the energy loss per half orbit with increasing energy seen in the two figures. In the Markovian 
case (t = 0), the energy loss increases with increasing e. Two opposing effects contribute to this 
behavior. On the one hand, the rate of energy loss is directly proportional to the kinetic energy 
of the particle and, on average, the kinetic energy increases with increasing e. On the other hand, 
the energy loss in a half orbit increases with the time spent in this half orbit, and this time, t e , 
decreases with increasing e. These competing effects more or less (but not entirely) compensate 
in the Markovian limit, and the result is a mild increase of fi with increasing e. The situation 
becomes more complicated in the non-Markovian case, that is, when t > 0. Now the overwhelming 
factor in the e dependence of the energy loss per half orbit is the relation between r and t e . More 
specifically, it is the relation between ttt and t £ , as seen from Eq. (|34|). In the regime ttt « t £ 
one observes a behavior similar to the Markovian case, that is, in this regime \x increases with e (in 
Eq. ( p^D this has been approximated by the upper bound discussed earlier, that is, by the constant 
47/3; in any case, this regime is not visible in Fig. |7|). However, when ttt > t £ the memory term is 



not fully dissipative (cf. discussion surrounding Eq. (17)) and hence leads to a smaller energy loss 
per half orbit. Since t e decreases with increasing e, jj, correspondingly decreases as well. Thus, for 
finite r the energy loss per half orbit as a function of energy above the barrier is expected to exhibit 
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Figure 7: Energy loss per half orbit as a function of energy in the non-Markovian regime for 
correlation time r = 3.0 and dissipation parameter 7 = 0.01. 

a maximum. The behavior seen in Fig. [?] is that of fi to the right of this maximum. Note that r 
would have to be extremely small in order to see the initial rise in fi (which in this approximation 
would be a flat region). Note also that although both Fourier components are needed to reproduce 
the behavior of p(t), the behavior of //(e) is determined overwhelmingly by the principal (the first) 
Fourier component in Eq. (|30|). 

From the result Eq. (|3~4|) we now wish to extract a specific (constant) value of fi to use in 
the calculation of the transmission coefficient. This is more difficult in the non-Markovian case 
than in the Markovian case because the variation of \i with e is much more pronounced. However, 
one should remember that at low temperatures the population above the barrier is small and the 
principal contributors to the transmission coefficient are the particles near the barrier. Indeed, 
in the Markovian case we see from Fig. |B| that there are about five or six oscillations before the 
transmission coefficient settles down, indicating the participation of particles in a layer of width 
5/j, or so above the barrier. The variation of fi over a comparable range in Fig. [7| is actually rather 
small, that is, only a small portion of the range shown in the figure is actually important for the 
transmission coefficient. Thus, for instance, for r = 3 we find that the value of /t for those particles 
that half-orbit once before becoming trapped (found from Eq. (|34|) by setting e = /t and solving 
for (j,) is /1/7 = 0.3175. For 7 = 0.01 this corresponds to \x = 0.003175. The variation of // between 
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e = 0.003175 and, say, five times this value, is about 15%. Thus, if we take fj, = 0.003175 we might 
overestimate \x by an amount only somewhat greater than that of the Markovian calculation, and 
expect the result to be comparably satisfactory. 
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Figure 8: Transmission coefficient of an ensemble of particles at temperature ksT = 0.025, 
dissipation parameter 7 = 0.01, and dissipative memory time r = 3.0. Solid curve: simulation of 
4000 particles. Dashed curve: our theory, Eq. ( |23| ) with fi = 0.003175. 

In Fig. H we show the simulation results as well as our theoretical result for a typical set of 
parameters. The simulation curve is one of those shown in Fig. ||. As in the Markovian case, the 
theory captures the distinctive features mentioned in the Markovian case: the time of the first 
abrupt drop of the transmission coefficient, the frequencies and amplitudes of the oscillations, and 
the asymptotic value. The theoretical overestimation of n in the later oscillations is discernible (as 
it is in Fig. |5|), and the calculated equilibrium transmission coefficient is also perhaps a bit high 
(as it is in Fig. |5|). However, we understand the causes of these small effects and know how to 
correct the theory for them if needed. Also, the thermal corrections would bring the equilibrium 
transmission coefficient into even closer agreement with the simulation results, as in the Markovian 
case. Altogether, the theory is clearly very good in its prediction of the complex time-dependent 
and equilibrium behavior of the transmission coefficient in the non-Markovian regime for the case 
of an exponential memory kernel. 
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4 Conclusions 



In this paper we have analyzed the time-dependent transmission coefficient for the generalized 
Kramers problem, that is, for a particle in a bistable potential described by a generalized Langevin 
equation. The transmission coefficient is associated with the rate of transitions of the particle from 
one well to the other. In this presentation the dissipative memory kernel has been taken to be 
exponential in time. In the diffusion-limited regime we have confirmed the predictions of Kohen 
and Tannor || via numerical simulations, and have shown their predictions to be accurate both in 
the non-adiabatic regime and in the caging regime. In the energy-diffusion-limited regime we have 
generalized our previous theory for the Markovian limit [|l0| to the non-Markovian case and have 
shown good agreement with numerical simulations in this regime as well. One can therefore say 
that there is now a fairly complete understanding of the time-dependent transmission coefficient 
for the generalized Kramers problem with an exponential memory kernel in the limiting regimes of 
strong dissipation and weak dissipation. 

These remarks cover three of the four regimes represented by the simulations in Fig. ^, that 
is, the 7 = 2.0, 7 = 10, and 7 = 0.01 curves. The fourth curve, the one for 7 = 0.3, displays a 
distinct behavior that has not been discussed in this paper but that will be discussed elsewhere 
p0| . This represents an "intermediate" dissipation and is thus associated with a regime in which 
elements appropriate to the diffusion-limited regime and also to the energy-diffusion-limited regime 
in combination play an important role. 

If one were to apply (misadvisedly) the KT theory to the left of the turnover in Fig. [I], one 
would find an apparent settling of the transmission coefficient at the values that in reality represent 
only the short-time plateau of K(t) seen in Fig. |2| before the first abrupt drop. Thus, for 7 = 0.01 
the transmission coefficient would settle at nearly unity and for 7 = 0.3 it would settle at around 
0.96. The actual abrupt drops seen in Fig. |2| following these initial plateaus come about because of 
the particles that have not been trapped after one half orbit and that at times t £ cross the barrier 
to move over the left well. The diffusion-limited theory of course does not capture this recrossing 
and hence levels off. This leveling off has been referred to as an example of a "false plateau" that 
one must be careful of in interpreting simulation data. Perhaps more interesting "false plateau" 
manifestations occur as a result of memory kernels that are not exponential and that can lead to 
bottlenecks of various kinds in the relaxation process [|l^]. In general, one expects that different 
memory kernels might lead to behaviors not seen in the particular exponential model considered in 
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this paper. Our analysis of other memory kernels and the resultant exploration of interesting new 
dynamical manifestations of these differences will also be presented elsewhere ^Cj . 
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Figure 1: Transmission coefficient K s t vs dissipation parameter 7 for ksT = 0.025 and various 
values of the memory correlation time r. These results are obtained from direct simulation of 
Eq. (|l]) as described in Section 2.2. As is well known, the turnover moves to higher values of the 



friction parameter with increasing r because the effective dissipation at given 7 is weaker as r 
increases. Also, the maximum transmission coefficient decreases with increasing r. 



Figure 2: Numerical simulation results for the transmission coefficient n{t) vs time for memory 
correlation time r = 3.0 and temperature k^T = 0.025, and four values of the dissipation parameter 
7. Diffusion-limited non-adiabatic regime: 7 = 2; diffusion-limited caging regime: 7 = 10; energy- 
diffusion-limited regime: 7 = 0.01; intermediate regime: 7 = 0.3. 



Figure 3: Numerical simulation results and theoretical results for the transmission coefficient n(t) 
vs time for memory correlation time r = 3.0 and temperature k^T = 0.025, and two values of 
the dissipation parameter 7 in the diffusion-limited regime. Solid curves: simulation results for 
7 = 2.0 (diffusion-limited non-adiabatic regime) and 7 = 10.0 (diffusion-limited caging regime). 
The dashed curves in each case are calculated from the theory of Kohen and Tannor as given in 
Eq. 



Figure 4: Energy loss per half orbit as a function of energy in the Markovian limit (r — * 0) for 
two values of the dissipation parameter: 7 = 0.005 (dashed curve) and 7 = 0.01 (solid curve). 



Figure 5: Transmission coefficient of an ensemble of particles at temperature fcgT = 0.025 and 
dissipation parameter 7 = 0.005 in the Markovian limit. Solid curve: simulation of 4000 particles. 
Dashed curve: our theory, Eq. (23), with the bare value [i = 47/3 = 0.006667. Thick line that 
intersects the right vertical axis: value of the equilibrium transmission coefficient n s t obtained from 
Eq. (|26|) using the value of fi renormalized by thermal fluctuations, which for this temperature is 
fi(T) = 0.004907. 



Figure 6: Momentum as a function of time for a particle moving in the potential O with energy 
e = 0.003175 above the barrier. Solid curve: exact result. Dashed curve: approximation (30) with 

D- 



Figure 7: Energy loss per half orbit as a function of energy in the Markovian limit (r — * 0) for 
two values of the dissipation parameter: 7 = 0.005 (dashed curve) and 7 = 0.01 (solid curve). 



Figure 8: Transmission coefficient of an ensemble of particles at temperature ksT = 0.025, 
dissipation parameter 7 = 0.01, and dissipative memory time r = 3.0. Solid curve: simulation of 
4000 particles. Dashed curve: our theory, Eq. ( |23| ) with \i = 0.003175. 
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